ODS 1
1.1.1
library(readxl)
url <- "https://github.com/ChiaraZamoraM/ODS/raw/main/1.1.1.xlsx"
destfile <- "X1_1_1.xlsx"
curl::curl_download(url, destfile)
X1_1_1 <- read_excel(destfile, skip = 28)
New names:
* `` -> ...1
* `` -> ...2
* `` -> ...3
* `` -> ...4
* `` -> ...5
* ...
colnames(X1_1_1) = paste0("ODS1_1_1_",colnames(X1_1_1))
X1_1_1 = X1_1_1[,c(1,10:26)]
names(X1_1_1)
[1] "ODS1_1_1_...1" "ODS1_1_1_...10" "ODS1_1_1_2013" "ODS1_1_1_...12" "ODS1_1_1_2014"
[6] "ODS1_1_1_...14" "ODS1_1_1_2015" "ODS1_1_1_...16" "ODS1_1_1_2016" "ODS1_1_1_...18"
[11] "ODS1_1_1_2017" "ODS1_1_1_...20" "ODS1_1_1_2018" "ODS1_1_1_...22" "ODS1_1_1_2019"
[16] "ODS1_1_1_...24" "ODS1_1_1_2020" "ODS1_1_1_...26"
X1_1_1[,c(2:18)]=lapply(X1_1_1[,c(2:18)], as.numeric)
Warning in lapply(X1_1_1[, c(2:18)], as.numeric) :
NAs introduced by coercion
Warning in lapply(X1_1_1[, c(2:18)], as.numeric) :
NAs introduced by coercion
Warning in lapply(X1_1_1[, c(2:18)], as.numeric) :
NAs introduced by coercion
Warning in lapply(X1_1_1[, c(2:18)], as.numeric) :
NAs introduced by coercion
Warning in lapply(X1_1_1[, c(2:18)], as.numeric) :
NAs introduced by coercion
Warning in lapply(X1_1_1[, c(2:18)], as.numeric) :
NAs introduced by coercion
Warning in lapply(X1_1_1[, c(2:18)], as.numeric) :
NAs introduced by coercion
Warning in lapply(X1_1_1[, c(2:18)], as.numeric) :
NAs introduced by coercion
Warning in lapply(X1_1_1[, c(2:18)], as.numeric) :
NAs introduced by coercion
Warning in lapply(X1_1_1[, c(2:18)], as.numeric) :
NAs introduced by coercion
Warning in lapply(X1_1_1[, c(2:18)], as.numeric) :
NAs introduced by coercion
Warning in lapply(X1_1_1[, c(2:18)], as.numeric) :
NAs introduced by coercion
Warning in lapply(X1_1_1[, c(2:18)], as.numeric) :
NAs introduced by coercion
Warning in lapply(X1_1_1[, c(2:18)], as.numeric) :
NAs introduced by coercion
Warning in lapply(X1_1_1[, c(2:18)], as.numeric) :
NAs introduced by coercion
Warning in lapply(X1_1_1[, c(2:18)], as.numeric) :
NAs introduced by coercion
Warning in lapply(X1_1_1[, c(2:18)], as.numeric) :
NAs introduced by coercion
X1_1_1$ODS1_1_1_2013= (X1_1_1$ODS1_1_1_2013+ X1_1_1$ODS1_1_1_...12)/2
X1_1_1$ODS1_1_1_2014= (X1_1_1$ODS1_1_1_2014+ X1_1_1$ODS1_1_1_...14)/2
X1_1_1$ODS1_1_1_2015= (X1_1_1$ODS1_1_1_2015+ X1_1_1$ODS1_1_1_...16)/2
X1_1_1$ODS1_1_1_2016= (X1_1_1$ODS1_1_1_2016+ X1_1_1$ODS1_1_1_...18)/2
X1_1_1$ODS1_1_1_2017= (X1_1_1$ODS1_1_1_2017+ X1_1_1$ODS1_1_1_...20)/2
X1_1_1$ODS1_1_1_2018= (X1_1_1$ODS1_1_1_2018+ X1_1_1$ODS1_1_1_...22)/2
X1_1_1$ODS1_1_1_2019= (X1_1_1$ODS1_1_1_2019+ X1_1_1$ODS1_1_1_...24)/2
X1_1_1$ODS1_1_1_2020= (X1_1_1$ODS1_1_1_2020+ X1_1_1$ODS1_1_1_...26)/2
names(X1_1_1)[1]= "DEPARTAMENTO"
X1_1_1 <- data.frame(X1_1_1[,seq(1,18,2)])
X1_1_1$DEPARTAMENTO= gsub("Lima Metropolitana","LIMA",X1_1_1$DEPARTAMENTO)
X1_1_1$DEPARTAMENTO= gsub("Lima","LIMA PROVINCIAS",X1_1_1$DEPARTAMENTO)
1.2.1
url <- "https://github.com/ChiaraZamoraM/ODS/raw/main/1.2.1.xlsx"
destfile <- "X1_2_1.xlsx"
curl::curl_download(url, destfile)
X1_2_1 <- read_excel(destfile, skip = 28)
New names:
* `` -> ...1
* `` -> ...2
* `` -> ...3
* `` -> ...4
* `` -> ...6
* ...
colnames(X1_2_1) = paste0("ODS1_2_1_",colnames(X1_2_1))
X1_2_1 = X1_2_1[,c(1,16:32)]
names(X1_2_1)
[1] "ODS1_2_1_...1" "ODS1_2_1_...16" "ODS1_2_1_2013" "ODS1_2_1_...18" "ODS1_2_1_2014"
[6] "ODS1_2_1_...20" "ODS1_2_1_2015" "ODS1_2_1_...22" "ODS1_2_1_2016" "ODS1_2_1_...24"
[11] "ODS1_2_1_2017" "ODS1_2_1_...26" "ODS1_2_1_2018" "ODS1_2_1_...28" "ODS1_2_1_2019"
[16] "ODS1_2_1_...30" "ODS1_2_1_2020" "ODS1_2_1_...32"
X1_2_1[,c(2:18)]=lapply(X1_2_1[,c(2:18)], as.numeric)
Warning in lapply(X1_2_1[, c(2:18)], as.numeric) :
NAs introduced by coercion
Warning in lapply(X1_2_1[, c(2:18)], as.numeric) :
NAs introduced by coercion
Warning in lapply(X1_2_1[, c(2:18)], as.numeric) :
NAs introduced by coercion
Warning in lapply(X1_2_1[, c(2:18)], as.numeric) :
NAs introduced by coercion
Warning in lapply(X1_2_1[, c(2:18)], as.numeric) :
NAs introduced by coercion
Warning in lapply(X1_2_1[, c(2:18)], as.numeric) :
NAs introduced by coercion
Warning in lapply(X1_2_1[, c(2:18)], as.numeric) :
NAs introduced by coercion
Warning in lapply(X1_2_1[, c(2:18)], as.numeric) :
NAs introduced by coercion
Warning in lapply(X1_2_1[, c(2:18)], as.numeric) :
NAs introduced by coercion
Warning in lapply(X1_2_1[, c(2:18)], as.numeric) :
NAs introduced by coercion
Warning in lapply(X1_2_1[, c(2:18)], as.numeric) :
NAs introduced by coercion
Warning in lapply(X1_2_1[, c(2:18)], as.numeric) :
NAs introduced by coercion
Warning in lapply(X1_2_1[, c(2:18)], as.numeric) :
NAs introduced by coercion
Warning in lapply(X1_2_1[, c(2:18)], as.numeric) :
NAs introduced by coercion
Warning in lapply(X1_2_1[, c(2:18)], as.numeric) :
NAs introduced by coercion
Warning in lapply(X1_2_1[, c(2:18)], as.numeric) :
NAs introduced by coercion
Warning in lapply(X1_2_1[, c(2:18)], as.numeric) :
NAs introduced by coercion
X1_2_1$ODS1_2_1_2013= (X1_2_1$ODS1_2_1_2013+ X1_2_1$ODS1_2_1_...18)/2
X1_2_1$ODS1_2_1_2014= (X1_2_1$ODS1_2_1_2014+ X1_2_1$ODS1_2_1_...20)/2
X1_2_1$ODS1_2_1_2015= (X1_2_1$ODS1_2_1_2015+ X1_2_1$ODS1_2_1_...22)/2
X1_2_1$ODS1_2_1_2016= (X1_2_1$ODS1_2_1_2016+ X1_2_1$ODS1_2_1_...24)/2
X1_2_1$ODS1_2_1_2017= (X1_2_1$ODS1_2_1_2017+ X1_2_1$ODS1_2_1_...26)/2
X1_2_1$ODS1_2_1_2018= (X1_2_1$ODS1_2_1_2018+ X1_2_1$ODS1_2_1_...28)/2
X1_2_1$ODS1_2_1_2019= (X1_2_1$ODS1_2_1_2019+ X1_2_1$ODS1_2_1_...30)/2
X1_2_1$ODS1_2_1_2020= (X1_2_1$ODS1_2_1_2020+ X1_2_1$ODS1_2_1_...32)/2
names(X1_2_1)[1]= "DEPARTAMENTO"
X1_2_1 <- data.frame(X1_2_1[,seq(1,18,2)])
X1_2_1$DEPARTAMENTO= gsub("Lima Metropolitana","LIMA",X1_2_1$DEPARTAMENTO)
X1_2_1$DEPARTAMENTO= gsub("Lima","LIMA PROVINCIAS",X1_2_1$DEPARTAMENTO)
1.3.1
url <- "https://github.com/ChiaraZamoraM/ODS/raw/main/1.3.1.xlsx"
destfile <- "X1_3_1.xlsx"
curl::curl_download(url, destfile)
X1_3_1 <- read_excel(destfile, skip = 7)
New names:
* `` -> ...1
colnames(X1_3_1) = paste0("ODS1_3_1_",colnames(X1_3_1))
X1_3_1[,c(2:14)]=lapply(X1_3_1[,c(2:14)], as.numeric)
names(X1_3_1)[1]= "DEPARTAMENTO"
X1_3_1$DEPARTAMENTO= gsub("Prov. Const. del Callao","Callao",X1_3_1$DEPARTAMENTO)
X1_3_1$DEPARTAMENTO= gsub("Provincia de Lima","LIMA",X1_3_1$DEPARTAMENTO)
X1_3_1$DEPARTAMENTO= gsub("Región Lima","LIMA PROVINCIAS",X1_3_1$DEPARTAMENTO)
1.4.1
ODS_1 = Reduce(function(x, y) merge(x, y, by= "DEPARTAMENTO"), list(X1_1_1, X1_2_1, X1_3_1, X1_4_1))
X1_4_1[,c(2:8)]=lapply(X1_4_1[,c(2:8)], as.numeric)
names(X1_4_1)[1]= "DEPARTAMENTO"
X1_4_1$DEPARTAMENTO= gsub("Prov. Const. del Callao","Callao",X1_4_1$DEPARTAMENTO)
X1_4_1$DEPARTAMENTO= gsub("Provincia de Lima","LIMA",X1_4_1$DEPARTAMENTO)
X1_4_1$DEPARTAMENTO= gsub("Región Lima","LIMA PROVINCIAS",X1_4_1$DEPARTAMENTO)
Merge
ODS_1 = Reduce(function(x, y) merge(x, y, by= "DEPARTAMENTO"), list(X1_1_1, X1_2_1, X1_3_1, X1_4_1))
ODS 2
url <- "https://www.inei.gob.pe/media/MenuRecursivo/publicaciones_digitales/Est/Lib1795/cuadros/Cap010.xls"
destfile <- "Cap010.xls"
curl::curl_download(url, destfile)
Cap010 <- read_excel(destfile,sheet = "10.18",skip = 6)
New names:
* `` -> ...2
* `` -> ...3
* `` -> ...4
* `` -> ...5
* `` -> ...6
* ...
names(Cap010)[c(1,4)]= c("DEPARTAMENTO","2020")
Cap010 =Cap010[,c(1,4)]
Cap010$DEPARTAMENTO =gsub(' [0-9].', '', Cap010$DEPARTAMENTO)
url <- "https://www.inei.gob.pe/media/MenuRecursivo/publicaciones_digitales/Est/Lib1795/cuadros/Cap010.xls"
destfile <- "Cap010.xls"
curl::curl_download(url, destfile)
Cap010_2 <- read_excel(destfile,sheet = "10.18",skip = 6)
New names:
* `` -> ...2
* `` -> ...3
* `` -> ...4
* `` -> ...5
* `` -> ...6
* ...
names(Cap010_2)[c(1,7)]= c("DEPARTAMENTO","2020")
Cap010_2 =Cap010_2[,c(1,7)]
Cap010_2$DEPARTAMENTO =gsub(' [0-9].', '', Cap010_2$DEPARTAMENTO)
2.1.1
url <- "https://github.com/ChiaraZamoraM/ODS/raw/main/2.2.1.xlsx"
destfile <- "X2_2_1.xlsx"
curl::curl_download(url, destfile)
X2_2_1 <- read_excel(destfile, skip = 7)
New names:
* `` -> ...1
names(X2_2_1)[1]= "DEPARTAMENTO"
Cap010$DEPARTAMENTO= gsub("Prov. Const. del Callao","Callao",Cap010$DEPARTAMENTO)
Cap010$DEPARTAMENTO= gsub("Lima Metropolitana","LIMA",Cap010$DEPARTAMENTO)
Cap010$DEPARTAMENTO= gsub("Departamento de Lima","LIMA PROVINCIAS",Cap010$DEPARTAMENTO)
X2_2_1$DEPARTAMENTO= gsub("Prov. Const. del Callao","Callao",X2_2_1$DEPARTAMENTO)
X2_2_1$DEPARTAMENTO= gsub("Provincia de Lima","LIMA",X2_2_1$DEPARTAMENTO)
X2_2_1$DEPARTAMENTO= gsub("Región Lima","LIMA PROVINCIAS",X2_2_1$DEPARTAMENTO)
X2_2_1 = merge(X2_2_1, Cap010, by = "DEPARTAMENTO")
names(X2_2_1)[c(2:15)] = paste0("ODS2_2_1_",names(X2_2_1)[c(2:15)])
2.2.1
url <- "https://github.com/ChiaraZamoraM/ODS/raw/main/2.2.2.xlsx"
destfile <- "X2_2_2.xlsx"
curl::curl_download(url, destfile)
X2_2_2 <- read_excel(destfile, skip = 7)
New names:
* `` -> ...1
names(X2_2_2)[1]= "DEPARTAMENTO"
Cap010_2$DEPARTAMENTO= gsub("Prov. Const. del Callao","Callao",Cap010_2$DEPARTAMENTO)
Cap010_2$DEPARTAMENTO= gsub("Lima Metropolitana","LIMA",Cap010_2$DEPARTAMENTO)
Cap010_2$DEPARTAMENTO= gsub("Departamento de Lima","LIMA PROVINCIAS",Cap010_2$DEPARTAMENTO)
X2_2_2$DEPARTAMENTO= gsub("Prov. Const. del Callao","Callao",X2_2_2$DEPARTAMENTO)
X2_2_2$DEPARTAMENTO= gsub("Provincia de Lima","LIMA",X2_2_2$DEPARTAMENTO)
X2_2_2$DEPARTAMENTO= gsub("Región Lima","LIMA PROVINCIAS",X2_2_2$DEPARTAMENTO)
X2_2_2 = merge(X2_2_2, Cap010_2, by = "DEPARTAMENTO")
colnames(X2_2_2)[c(2:12)] = paste0("ODS2_2_2_",colnames(X2_2_2)[c(2:12)])
str(X2_2_2)
'data.frame': 34 obs. of 12 variables:
$ DEPARTAMENTO : chr "Amazonas" "Áncash" "Apurímac" "Arequipa" ...
$ ODS2_2_2_2010: chr "0.6" "0.3" "0.3" "0.5" ...
$ ODS2_2_2_2011: chr "0.8" "0.3" "0.5" "0.3" ...
$ ODS2_2_2_2012: chr "0.7" "0.4" "0.5" "0" ...
$ ODS2_2_2_2013: chr "0.1" "0.3" "0.4" "0.8" ...
$ ODS2_2_2_2014: num 0.4 0.2 0.3 0.3 0.2 0.4 0 0.3 0.2 0.9 ...
$ ODS2_2_2_2015: num 0.7 0.6 1.2 0 0.9 0.6 0.5 0.4 1.8 1.4 ...
$ ODS2_2_2_2016: num 0.6 0.1 0.7 0 1.2 0.4 0.2 1.7 1.3 0.3 ...
$ ODS2_2_2_2017: num 1.3 0.1 0.6 0.7 0.7 0.5 0.2 0.8 0.3 0.6 ...
$ ODS2_2_2_2018: num 1.2 0.2 0.4 0.1 0.5 0.2 0 0.3 0.7 0.3 ...
$ ODS2_2_2_2019: num 1.1 0.1 0.1 0.1 0.6 1 0.2 0.8 0.2 0.1 ...
$ ODS2_2_2_2020: num 0.9385 0.0757 0.3944 0.1938 0.6927 ...
Merge
ODS_2 = Reduce(function(x, y) merge(x, y, by= "DEPARTAMENTO"), list(X2_2_1,X2_2_2))
ODS_2[,c(2:26)]=lapply(ODS_2[,c(2:26)], as.numeric)
Warning in lapply(ODS_2[, c(2:26)], as.numeric) :
NAs introduced by coercion
Warning in lapply(ODS_2[, c(2:26)], as.numeric) :
NAs introduced by coercion
Warning in lapply(ODS_2[, c(2:26)], as.numeric) :
NAs introduced by coercion
Warning in lapply(ODS_2[, c(2:26)], as.numeric) :
NAs introduced by coercion
Warning in lapply(ODS_2[, c(2:26)], as.numeric) :
NAs introduced by coercion
Warning in lapply(ODS_2[, c(2:26)], as.numeric) :
NAs introduced by coercion
Warning in lapply(ODS_2[, c(2:26)], as.numeric) :
NAs introduced by coercion
Warning in lapply(ODS_2[, c(2:26)], as.numeric) :
NAs introduced by coercion
Warning in lapply(ODS_2[, c(2:26)], as.numeric) :
NAs introduced by coercion
Warning in lapply(ODS_2[, c(2:26)], as.numeric) :
NAs introduced by coercion
Warning in lapply(ODS_2[, c(2:26)], as.numeric) :
NAs introduced by coercion
Merge
ODS = Reduce(function(x, y) merge(x, y, by= "DEPARTAMENTO"), list(ODS_1,ODS_2))
library(stringi)
ODS$DEPARTAMENTO= stri_trans_general(str = toupper(ODS$DEPARTAMENTO), id = "Latin-ASCII")
ODSTrans$Ano = year(ODSTrans$Ano)
Error in as.POSIXlt.character(x, tz = tz(x)) :
character string is not in a standard unambiguous format
mapa_prov <- st_read("PROVINCIAS.shp")
Reading layer `PROVINCIAS' from data source `C:\Users\soyma\Documents\GitHub\ODS\PROVINCIAS.shp' using driver `ESRI Shapefile'
Simple feature collection with 196 features and 6 fields
Geometry type: MULTIPOLYGON
Dimension: XY
Bounding box: xmin: -81.32823 ymin: -18.35093 xmax: -68.65228 ymax: -0.03860597
Geodetic CRS: WGS 84
summary(subset1_1_1$Valor)
Min. 1st Qu. Median Mean 3rd Qu. Max.
0.25 0.50 3.05 4.29 6.85 20.25
subset1_2_1= subset(mapaODS, ODSNro=="1_2_1" & Ano> 2014)
base1_2_1= ggplot(data = subset1_2_1) + theme_light()
library(plotly)
mapaley2 = base1_2_1 +
geom_sf(aes(fill= Valor)) + labs(fill = "Porcentaje (%)") +
geom_sf(data = mapa_prov,
fill = NA) +
facet_wrap(~Ano)
mapa2 = mapaley2 +
scale_fill_gradient(low = "yellow", high = "darkred")+
labs(title = "Incidencia de la pobreza extrema", subtitle = "Indicador 1.1.1")
mapa2
subset1_3_1= subset(mapaODS, ODSNro=="1_3_1" & Ano> 2014)
base1_3_1= ggplot(data = subset1_3_1) + theme_light()
library(plotly)
mapaley3 = base1_3_1 +
geom_sf(aes(fill= Valor)) + labs(fill = "Porcentaje (%)") +
geom_sf(data = mapa_prov,
fill = NA) +
facet_wrap(~Ano)
mapa3 = mapaley3 +
scale_fill_gradient(low = "lightpink", high = "firebrick2")+
labs(title = "Proporción de población de 14 a más años de edad \ncon seguro de pensión", subtitle = "Indicador 1.3.1")
mapa3
subset1_4_1= subset(mapaODS, ODSNro=="1_4_1" & Ano> 2014)
base1_4_1= ggplot(data = subset1_4_1) + theme_light()
library(plotly)
mapaley4 = base1_4_1 +
geom_sf(aes(fill= Valor)) + labs(fill = "Porcentaje (%)") +
geom_sf(data = mapa_prov,
fill = NA) +
facet_wrap(~Ano)
mapa4 = mapaley4 +
scale_fill_gradient(low = "lightpink", high = "firebrick2")+
labs(title = "Proporción de la población que vive en hogares con \nacceso a servicios básicos de infraestructura", subtitle = "Indicador 1.4.1")
mapa4
library(leaflet)
pal1 =colorNumeric(palette = "OrRd", domain = subset1_1_1$Valor])
Error: unexpected ']' in "pal1 =colorNumeric(palette = "OrRd", domain = subset1_1_1$Valor]"
names(subset1_1_1)
[1] "DEPARTAMEN" "IDDPTO" "id" "ODS_Ano" "Valor" "Ano" "ODSNro"
[8] "geometry"